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We consider the temperature dependence of the sound attenuation and sound velocity in layered 
impure metals due to s-wave superconducting fluctuations of the order parameter above the criti- 
cal temperature. We obtain the dependence on material properties of these fluctuation corrections 
in the hydrodynamic limit, where the electron mean free path I is much smaller than the wave- 
length of sound and where the electron collision rate r~ is much larger than the sound frequency. 
For longitudinal sound propagating perpendicular to the layers, the open Fermi surface condition 
leads to a suppression of the divergent contributions to leading order, in contrast with the case of 
paraconductivity. The leading temperature dependent corrections, given by the Aslamazov-Larkin, 
Maki-Thompson and density of states terms, remain finite as T — > T c . Nevertheless, the sensitivity 
of new ultrasonic experiments on layered organic conductors should make these fluctuations effects 
measurable. 

PACS numbers: 74.25.Ld, 74.25.Fy, 74.40.+k 



The problem of fluctuations of the superconducting or- 
der parameter above the critical temperature began to at- 
tract the attention of researchers more than three decades 
a g°- CI- Q Considerable work has been done from both 
the theoretical and experimental points of view. Conduc- 
tivity, thermoconductivity, magnetoconductivity as well 
as tunneling properties and nuclear magnetic resonance 
characteristics are amongst the various transport phe- 
nomena where fluctuation contributions were predicted 
and experimentally observed. Ref. contains an exten- 
sive review of the field and references. 

In the present paper, we concentrate on the super- 
conducting fluctuation corrections to the sound veloc- 
ity and sound attenuation in layered superconductors 
above T c . As usual, superconducting fluctuations man- 
ifest themselves in three different ways. a) The ef- 
fective number of normal carriers is reduced because 
some of the electrons exist as transient Cooper pairs. 
This is the so-called density of states contribution, b) 
The single-particle excitations are Andreev reflected off 
the superconducting fluctuations, as described by the 
so-called Maki-Thompson term, c) Some of the elec- 
trons behave like Cooper pairs for a time given by the 
Ginzburg-Landau time. This is the famous Aslamazov- 
Larkin contribution. Layered organic superconductors, 
in particular, are an ideal class of materials to study 
fluctuation phenomena. This is because the relatively 
low charge-carrier concentration and the strong elec- 
tronic anisotropy enhance the effect of superconduct- 
ing fluctuations and increase the size of the fluctuation 
regime. The manifestation of superconducting fluctu- 
ations in organic materials was found in a number of 
experimental works on dc-magnetoconductivity, thermo- 
conductivity, ac-susceptibility, specific heat and torque 
magnetometry. || Preliminary ultrasonic experiments 
in the K-phase layered organic materials |(| (BEDT- 



TTF) 2 Cu[N(CN) 2 ]Br, (BEDT-TTF) 2 Cu[N(CN) 2 ]Cl and 
(BEDT-TTF) 2 Cu(NCS) 2 suggest that the sound attenu- 
ation coefficient is decreased by superconducting fluctu- 
ations over a sizeable temperature range above T c . How- 
ever, the effect is much weaker than the fluctuation cor- 
rection to the conductivity Q where, in the same com- 
pounds, the Aslamazov-Larkin contribution dominates 
the superconducting fluctuation effects. Although sound 
attenuation experiments seem to be a very promising 
method to study fluctuation phenomena in layered super- 
conductors, the theoretical description is still incomplete 
and needs to be analyzed in more details. This more de- 
tailed analysis can tell us, for example, which material 
properties control fluctuation phenomena. Fluctuation 
phenomena can then be used to access these properties 
when they are not independently known. For example, 
we will see that the phase breaking time can be very 
important to determine the fluctuation contribution to 
ultrasonic attenuation whereas it would not contribute 
to dominant terms in the corresponding fluctuation con- 
ductivity. It is by considering a large set of experimental 
probes that one can make sure that all material proper- 
ties are unambiguously determined. 

In this paper we assume, as usual for fluctuation cor- 
rections, that we are close enough to T c that the char- 
acteristic pair frequency is less than temperature (in di- 
mcnsionless units). This is the so-called renormalized- 
classical regime. We study the case where longitudinal 
sound propagates perpendicular to the layers. This is 
the most relevant case to study since it is most easily 
accessible experimentally. Q We also consider the hy- 
drodynamic limit where the electron mean free path I is 
much smaller than the wavelength of sound and where 
the electron collision rate r _1 is much larger than the 
sound frequency. In quasi two-dimensional systems, the 
Fermi surface is opened in the direction perpendicular 
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to the layers. The coupling to the sound comes from 
the modulation of the hopping integral across the layers, 
taking also into account the electroneutrality condition in 
the moving frame. Q The Frohlich model where phonons 
couple to the electronic density is not valid. The most 
dramatic consequence of the open Fermi surface in such 
a situation is the suppression of the divergence in the 
Aslamazov-Larkin and Maki-Thompson terms that of- 
ten give, by contrast, the most important contribution 
to the fluctuation correction to the electrical conductiv- 
ity. Hence, sound propagation is less strongly modified 
by superconducting fluctuations than the electrical con- 
ductivity. All the calculations are performed for s— wave 
superconducting fluctuations. The d— wave case needs a 
separate treatment. 0, 

In the following section, we discuss earlier work and 
present the Hamiltonian. The diagrams that enter the 
calculation and then details of the calculation appear in 
the following sections. The results are gathered in SccM 
and discussed in Sec lVIl The appendix contains mathe- 
matical details on the analytic continuation of digamma 
functions. 



I. THE MODEL AND DEFINITIONS 

Ultrasonic attenuation in impure metals was analyzed 
by Pippard [Tlj . using the Boltzmann equation method, 
and by Tsuneto Il2j with the density matrix formalism. 
Later on, Schmid |l3j . using Green function methods, 
derived the effective electron-phonon interaction in the 
presence of impurities taking into account the screening 
of Coulomb repulsion between electrons. A key point of 
the approach of Tsuneto, Pippard and Schmid was to 
consider the electron system in a reference frame moving 
together with the ion lattice. In an impure metal, the 
elastic scattering of electrons as well as perfect screening 
at small phonon momentum and relaxation to equilib- 
rium occur in this oscillating frame. In this formalism, 
the electron-phonon coupling appears through the stress 
tensor instead of through the density operator occurring 
in the Frohlich model. The latter model leads to erro- 
neous results in the dirty limit. 0] 

The idea of a moving reference frame was also used 
by Kotliar and Ramakrishnan 0|, who considered the 
limit of strongly disordered metal and analyzed the effect 
on ultrasound propagation of incipient localization near 
the metal- insulator transition. Reizer considered the 
effect of various types of inelastic scattering (including 
electron-electron, electron-magnon scattering and weak- 
localization effects) on sound attenuation. 

The fluctuation corrections to the ultrasound attenu- 
ation in 3D metals were estimated in the most singular 
channel in the early work of Aslamazov and Larkin. [J 
Aslamazov and Varlamov, back in 1979, found the fluc- 
tuation corrections in layered superconductors. How- 
ever, these papers did not consider explicitly the form of 
the density of states for a corrugated cylindrical Fermi 



surface. That density of states is energy independent 
contrary to the assumption made in these papers. In ad- 
dition, the calculation was done for the Frohlich model. 
Also, Ref. 0| used impurity vertex corrections for the 
external vertices in the correlators and we shall see that 
these are unnecessary in the hydrodynamic limit when 
one takes into account the fact that the calculation of 
impurity scattering should be done in the electrically neu- 
tral moving frame. [l3l IT3| For the case of interest, with 
open Fermi surface in the perpendicular direction, we find 
results that are quite different from these earlier works. 
We use the unperturbed quasiparticle energy spectrum 
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C(P||.P*) = -\~ h 2t_LCOs(p z a), (1) 

appropriate for layered metals. Here, p z is the component 
of the quasiparticle momentum perpendicular to the con- 
ducting planes, pii is the in-plane momentum, t± <C Ef 
is the small intcrlaycr hopping energy, and a is the in- 
tcrlaycr distance. In this model, the Fermi surface has 
the well-known shape of a corrugated cylinder, and the 
density of states is given by uq — m/(2ira). It is impor- 
tant to note that in this model the density of states Vq 
does not depend on the quasiparticle energy. In other 
words, the derivative (di/o/d^) = vanishes. This rela- 
tion is obviously true at t± = 0, i.e. for non-interacting 
two-dimensional layers. If there is a small intcrlaycr in- 
teraction ij_ 7^ 0, any cross-section of the Fermi surface 
made at any value of p z still gives the two-dimensional 
(2D) result with (dv^/d^) = 0. One can easily see that 
this relation is valid for any Fermi surface that is opened 
in the third direction (t± <C £f) an d has circular cross 
section in the plane. 

As shown by Tsuneto |l2ll and Schmid 0], based on 
earlier work of Pippard |ll| and Blount, ^3] the canonical 
transformation that takes us to the moving frame leads, 
in the continuum limit, to an electron-phonon interaction 
that originates from the commutator of the infinitesimal 
generator of the transformation with the kinetic energy 
operator. The remaining contributions from the com- 
mutator of the generator of the canonical transformation 
with the electron-ion and electron-electron interactions 
are overall negligible. Since a trace over fermions is per- 
formed to compute phonon damping, the result should be 
independent of the canonical transformation. The mov- 
ing frame is the one best suited for approximations. 

For a tight binding model, such as the one we need 
to consider to obtain the proper open Fermi surface for 
propagation in the perpendicular direction p z , Eq. 
the corresponding microscopic derivation of the stress 
tensor has not been performed yet. It is not our pur- 
pose to give this derivation here. Instead, we follow the 
footsteps of Barisic, Su-Schrieffer-Heeger .JT9| and, 
more recently, Walker, Smith, and Samokhin [8j and as- 
sume that the electron-phonon interaction comes from 
the modulation of the hopping parameter t± induced by 
the lattice deformation. Neglecting umklapp processes, 
one obtains for the interaction of a longitudinal phonon 
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propagating in the z direction perpendicular to the planes 



(Huj fk) \ 
N ° M 2 ) (k' 2 ) ( C0S P* a ) c p+k,<T C p^ (al k + ak) , (2) 
s J 



where loq (k) = v s k is the sound frequency, v s the sound 
velocity, M the ion mass, AT the number of unit cells, 
G a constant that depends on the derivative of the hop- 
ping integral t± with respect to the strain, while a k ^ , Cp^ 
are, respectively, destruction and creation operators for 
phonons and for electrons of spin a. The above expres- 
sion is restricted to longitudinal phonons propagating in 
the perpendicular direction, hence the wave vector and 
polarization of the phonons involved in the interaction 
will satisfy kz = ez = 1. The above expression neglects 
the compression and stretching of chemical bonds that 
are not strictly along the direction of propagation. 

II. SOUND ATTENUATION WITHOUT 
FLUCTUATION CORRECTION 

The sound attenuation coefficient is determined by the 
imaginary part j(k) of the complex frequency u)(k) where 
the pole of the phonon Green function £)(k, w) is located. 
This quantity obeys Dyson's equation |20j 

D-\k,u„) = (zAk^r 1 -n(k,w„). (3) 

Expressed in bosonic Matsubara frequencies lo v = 2irvT 
using units fog = 1 H = 1, with v an integer, the quantity 
D°(k,w„) 

g°(k,*„) = - ffl (4) 

is the phonon propagator in the non-interacting case and 
IT(k,w t/ ) is the phonon self-energy. In an impure metal, 
the diagrammatic representation of n(k, ui v ) is given by 
Fig. 2] In the theory of electronic conductivity, one 
would consider a similar diagram for the current-current 
correlator with the bold dots representing the compo- 
nents of the vector velocity ev a = e{d^/ dp a ), where £ 
is the electronic energy and p a the a component of mo- 
mentum. In the case of sound attenuation, the vertex 
comes instead from the stress tensor corresponding to 
the electron-phonon interaction Eq. Using Eqs. @ 
and (@J), one obtains for this vertex gcos(p z a) where g is 
a constant. 

Each solid line in Fig. ^ corresponds to the electron 
Green function which, at finite temperature in the pres- 
ence of impurities, is given by 

G(p ' £n) ^Ii^?(F)' (5) 



p+k, £ W +G>v 




P, e« 



FIG. 1: The phonon self-energy in the case of an impure 
metal. Solid lines are one-electron Green functions, dotted 
lines the bare phonon propagators, dashed lines with crosses 
correspond to the impurity scattering, and bold circles show 
the renormalized electron-phonon vertices. 



where £(p) is the quasiparticle energy and where we de- 
fined e n — e n + l/(2r)sign(e„) with e n = irT(2n + 1) the 
Matsubara frequency and r the elastic scattering time. 
The quantity i/(2r)sign(e n ) is the imaginary part of the 
quasiparticle self-energy. The real part of the self-energy 
is constant and is absorbed in the definition of the chem- 
ical potential. 

In principle, one should include the vertex corrections 
that are represented, to leading order in (pf^) _1 with 
£ the electron mean free path, by the non-intersecting 
dashed impurity lines, as shown in Fig. ^ I n the case of 
the current-current correlator, vertex corrections lead to 
the replacement of the scattering time r by its transport 
analog r tr . [2(| For s-wave scattering, vertex corrections 
vanish because the vector vertex ev a averages to zero 
upon angular integration at vanishing external momen- 
tum k. This argument generally does not work with the 
correlator of the stress tensor. However, taking into ac- 
count perfect screening, it was shown by Schmid [l3T ] in 
the continuum limit that there is no impurity diffusion 
enhancement of the electron-phonon vertex in the case 
of transverse phonons, and that for longitudinal phonons 
this effect is negligible in the hydrodynamic limit k£ <C 1 
and wr<l. Physically, this comes about from the fact 
that the calculation should be done in the moving frame 
and that screening is perfect at long wavelengths. In our 
case, the analogous argument of electroncutrality leads 
@ to the replacement of the stress vertex F by F — (F) , 
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where the average accounts for the chemical potential 
shift. That average is precisely what is needed to 
make the impurity vertex correction vanish. In addition, 
in our specific case, (F) = {g cos(p z a)) = to the order 
in phonon wave vector k that we need. 

The phonon self-energy without fluctuation correction 
may finally be obtained from 



n(fc,w v ) 



(0) 



d 3 p 

(2^)3 



cos 2 (p z a)G(s n ,p)G(e n +w„,p + k). 



To evaluate this integral when r is finite, it is easiest to 

2 _ 2 

make the change of variables £ = -\ — - + 2tj_ cos(p z a) 



d 3 p m 
(2tt) 3 = (2tt) 3 



2,t 



(10 



7r / a 



7r / a 



d Pz (7) 



and to evaluate the integral over £ first, being careful to 
add and subtract the clean limit result to insure conver- 
gence. After summation over e n one should make the 
analytic continuation of the external phonon frequency 
following the rule iio v uj + iS. For sound, it is justi- 
fied to work in the hydrodynamic limit lot <§C 1 even for 
very pure systems. The quantity kvpT = kl also appears 
in the loop integral. Since the Fermi velocity is usually 
much larger than the sound velocity, kvp 3> lo — v s k, the 
limit kt > 1 may often be reached in ultrasonic experi- 
ments even though we always have lot -C 1. Hence, 
expansion of the result in powers of ki may not be justi- 
fied. In our case, however, since we are considering sound 
propagation perpendicular to the layers, vf is reduced 
by the anisotropy ratio t±/tu, which can be as small as 
10~ 4 , so that it is justified in our case to consider the 
limit fc!< 1. One can thus expand in powers of lot and 
kl. The terms with odd powers of lo give corrections to 
the sound attenuation while the terms proportional to 
even powers of lo contribute to the sound velocity. The 
first corrections that involve k appear to order (lot) (k£) 
because of the p z integral. In all that follows then, we 
set k = and expand in powers of (lot). 

—g 2 VQLUT one can obtain the con- 
tribution to the imaginary part of the phonon frequency 

M 
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7 H = -Loo(k) Im (U r (lo)) , 



(8) 



where is located the pole of Dyson's equation ®. We 
assume that the real part of the frequency at which the 
pole is located, lo, is close to the unperturbed frequency 
so that u>o(k) ~ u>o(k) 2 /lo (k). The power attenuation 
then reads 



a(uj) = —2j((jj)/v s 



(9) 



The renormalization of the phonon frequency Lo(k) 
is obtained from the real part of phonon self-energy 

Re [U(u, k)] R = -g 2 v a (l - (w) 2 ) using 



L0(k) = LO (k)Vl+RcU R - 



(10) 



III. DIAGRAMS FOR FLUCTUATION 
CONTRIBUTION TO SOUND ATTENUATION 

The preceding discussion shows that one can consider 
the same set of dia gram s, illustrated in Fig. |2J as for 
the conductivity pj. |22^. substituting gcos(p z a) for the 
vector vertices ev a at the extreme left and right-hand 
sides. Here, each wavy line in the figure corresponds 



1 js fs ~ r ^ r t>^ 










FIG. 2: Feynman diagrams giving the leading-order correc- 
tions from the fluctuations to the sound attenuation. Diagram 
1 is of the Aslamazov-Larkin (AL) type, diagrams 2-4 are of 
Maki- Thompson (MT) type, and the diagrams 5-10 are of 
the density of states (DOS) type. Solid lines are the normal 
state Green functions, wavy lines the fluctuation propaga- 
tors, shaded semicircles and shaded rectangles the impurity 
ladder averagings, dashed lines with cross the single impurity 
scatterings, and bold circles the renormalized electron-phonon 
vertices. Adapted from Ref.Q. 

to the fluctuation propagator (Cooper ladder) L(q, ilk) 
which, as T — ► T c , has the form: 



L(q,n k 



((A£(g,p) b | =p J 2 ) f .s. 
(47rT) 2 



in, 



4ttT J 



(11) 
-1 



where v s is the sound velocity. The sound attenuation 
increases with r, as does conductivity. 



where e f=a (T—T c )/T c , and tp(z) is the digamma function. 
In the above Eq. I|ll|) , the angular average over the Fermi 
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surface with the spectrum in Eq. Q reads: 
<(A£(q,p) b | =pp ) 2 >F.s. 



(12) 



= \ ((^9||) 2 + (4t± sin fea/2)) 2 ) = r^S*/ 2 , 

where D is a generalized diffusion operator. 

Returning to Fig. [3 the shaded semicircles correspond 
to vertex corrections from impurity averaging and are 
given by: 



A(q,£i,e 2 ) 



El - E 2 



£2\ 



Dq 2 



T* £l-£2 



r e(-£ie 2 ) 



(13) 



Diagrams 3, 4, 9 and 10 of Fig. [21 contain the impurity 
ladder in the particle-particle channel (shaded rectangle) : 



r c (q;£i,£ 2 ) 



2tttvq 



(14) 



1 (l 

£i -Ear V 



Dq 2 



r(£i - £ 2 ) 5 



while diagrams 7 and 8 have a single impurity line, which 
corresponds to the factor \/(2-ktuq). One can see that all 
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other insertions of impurity lines cither lead to diagrams 
that vanish, or give negligible corrections. 



IV. CALCULATIONAL DETAILS FOR THE 
DIAGRAMS 



A. MT and DOS diagrams 



We illustrate the procedure to evaluate the diagrams 
in Fig. by finding the expression for diagram 2 (Maki- 
Thompson type). In analytical form it is given by the 
integral: 



(15) 



2 -9 2t E / 4^L(q t Q k )K(q,!l kt u v ), 



where 



K{q,n k ,u v ) = ry^A(g,£n,O fc - e n )X(q,e n + u> v ,Ct k ~ e„ - u„) 
d 3 p 



(16) 



(2tt) 3 



cos (p z a) cos ((q z - p z )a) G(p, e n )G(p + k, e n + u) v )G{q -p,Q k - e n )G(q - p - k,fl k - e n - co v ), 

I 



For simplicity, and as in experiment, jo] we assume that 
the sound propagates in the direction perpendicular to 
the conducting layers, so that k has only one compo- 
nent: k = (0, k z ). The most singular contribution from 
the fluctuation propagator L(q, fl k ) comes from q — > 0. 
Therefore, one can neglect the q-dependence in the elec- 
tron Green functions. We also set fl k = since super- 
conducting fluctuations arc important in the renormal- 
ized classical regime where the characteristic frequency 
for the pair propagator is less than temperature. 

Performing the integration over p in Eq. H16J) first, one 
has for K(q, u> v ) = K(q, 0, u> u ): 

K(q, oj v ) — mcos(q z a)T (17) 
1 



E 



i 



i 



2|e, 



Dq 2 ) 2|£„|+ZV |£ 



2ty: 



1 



1 



1 



^2£„ 



Dq 2 2e n + Dq- 



2 \ 2e T 



T 



1 



- v \ 2e 



2u„ 



Dq 2 ) (- 



2£„ + Dq' 



Each of these two contributions leads to a different 
temperature dependence. The limits of summation in the 
second part of Eq. H18f) can be obtained by taking into 
account the fact that frequencies e„ = -kT(2u + 1) are of 
the fermionic type while the sound frequency lj v = 2-kTv 
is bosonic. The summation over e n in the second part of 
Eq. (jT%|) results in: 



m cos(q z a) 



4ttT I 



Dq 2 +2ul l 
AirT 



that is then divided into two parts. Namely, the first 
part has the summation in the interval e n G (-co, — lu u )U 
[0, oo), while in the second one the summation is done in 
the interval e n G [— uj v ,Q), 

K(q,cj„) (18) 



4tt(1 + w v t) 



(Dq 2 



(19) 



The analytical continuation of this sum involves some 
subtleties that are discussed in the Appendix. 

One can see that as T — > T c , the main g-dependence 
in Eq. i|15|) comes from the fluctuation propagator L(q) 
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and from the poles of K(q,u> v ). Therefore, it generally 
suffices to expand the regular part of K(q,Lu„) in powers 
of q and to keep the very first non- vanishing terms of the 
expansion. 

The singular part of Eq. (|19|l comes from taking the 
g°-term of the numerator together with the diffusion pole 
\/{Dq 2 + u;„) and analytically continuing them to real 
phonon frequencies. One arrives at integrals of the type: 



d 3 c 



cos(q z a) 



( 27r ) 3 (iV + 1/t*) (e + ml + r sin 2 (g z a/2)) ' 

(20) 

where r = 16t^_r]/vp <C 1 is the anisotropy parameter 
|23| . and 7* is the phase-breaking time, which provides 
the convergence of the integrals over the pair momentum 
q at small wave vector q. [3j The coefficient rj has the 
meaning of the square of the effective coherence length £ 
in the isotropic 2D case: 3 



v = e (tt) 



(tvfY 
2 



(21) 



Attx J Y \ 2 



1 , (\ 



At T = T c , the anisotropy parameter can be written as 
r = 4£ 2 _(0)/a 2 where £j_(0) is the Cooper pair size in the 
perpendicular (z) direction. We will also use the phase- 
breaking parameter 7^ = 2rj/(v 2 F tt*) instead of t*. 

The integral in Eq. (|20|l results in the so-called anoma- 
lous MT contribution. Note however that while the lead- 
ing term in the expansion of the numerator of Eq. 1(19(1 
gave the above result, Eq. (|2U|I . the q 2 -term of the expan- 
sion cancels the diffusion pole. This "regular anomalous" 
term, whose analog was overlooked in early work on con- 
ductivity fluctuation corrections, [24| results in the same 
type of integral as that given by the first part of Eq. ijT%|> 
with q = 0: 



d 3 c 



COs(q z a) 



(2n) 3 e + rjq 2 + r sin 2 (q z a/2) ' 



(22) 



Together, these two terms give so-called regular MT con- 
tribution, which generally has weaker temperature de- 
pendence than the anomalous MT part. 

The evaluation of the diagrams 5-8 of Fig. results in 
the integrals of the type 



(2tt)s 



rsin^(g z a/2)' 



(23) 



which give the density of states part of the fluctuation 
corrections. 

The remaining q- integration in Eqs. (|20|l . 1(22(1 . and 
(123(1 is similar to that encountered for the conductivity 
calculations and is described in much details in the re- 
view Q and in Ref. j^- We should emphasize that as 



T — > T c , the main temperature dependence of the fluctu- 
ation diagrams comes from these integrals. The remain- 
ing "bubble" , which consists of Green functions and im- 
purity lines, has a weaker dependence on Tt. It provides 
the coefficient in front of the function of e. 

In the conductivity calculation, one can drop the di- 
agrams containing the impurity ladder (3, 4, 9 and 10 
of Fig. I2J). This is because they contain the integration 
of one vector vertex with the Green function triangle. 
Angular integration over p leaves only terms containing 
powers of the small factor (v^fcr). In the general case of 
sound attenuation and propagation these diagrams would 
contribute, but for sound propagating perpendicular to 
the layers the integration over the open Fermi surface of 
the vertex gcos(p z a) also leads to a vanishing contribu- 
tion for diagrams 3, 4, 9 and 10 of Fig.[21to leading order 
in kt 

The small lu expansions of the DOS diagrams 5, 6, 
and of the regular contributions to the MT diagram 2 
each lead to zeroth-order terms, the combination of which 
does not sum to zero. One notices the difference with 
the conductivity calculation, where such a cancellation 
of zeroth-order terms ensures that there is no anomalous 
diamagnetism above T c . The DOS diagrams 7, 8, and 
the anomalous part of the MT diagram 2 do not have 
zeroth-order terms that would contribute to the leading 
renormalization of the sound velocity. 



B. Aslamazov-Larkin diagram 

Finally, let us consider the Aslamazov-Larkin (AL) di- 
agram 1, which gives the main fluctuation correction to 
most of the transport phenomena. In the case of sound 
propagation perpendicular to the layers, the analytical 
expression for the AL diagram reads: 

n^K,*) = -2.9 2 T£ / -0^B 2 (q,Q k ,u; v ) 



x L(q,n k )L(q,fl k +u v ), 



(24) 



where 



B(q,Q k ,Lo v ) = X(q,E n + w v ,—s n )X(q,e n ,-Sn)(25) 
d 3 p 



■ cos(p z a)G(p, e n )G{p, e n + w v )G(q - p, —e n ), 



In Eq. (|25|) . we have neglected the fife-dependence of 
the electron Green functions. One can also set w v = 
and expand G(q—p, —£„) in Eq. Ij25(l . [22j For sound prop- 
agating perpendicular to the planes, we can use r<Cl (See 
Eq .1)42(11 to find the appropriate expansion parameter for 
G(q — p, — e n ) appearing in the integral, namely: 

£(q P) £(P) 2 *i- (cos((g a - p z )a) - cos(p z a)) 



max{7rT, 1/r} 



max{7rT, 1 /r} 
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After the angular integration in Eq. (|25|l . the leading 
nonvanishing term reads: 



B(q) 



4tj_7]is a (l - cos(q z a)) 



(26) 



Following the Ref. one can integrate Ecr. (|24[l . which 
after the analytic continuation and expansion in powers 
of cj, results in the integral of the type: 



/ 



d 3 q 



sm 4 (q z a/2) 



(2?r)3 {e + m l+rsm\q z a/2j) 



n+2 



(27) 



in the n-th order of aj-expansion. To leading orders in 
oj, the integrals in Eqs. (|27|l arc not divergent at e — > 0. 
The resulting AL terms are finite and comparable with 
other fluctuation contributions. 



V. RESULTS 

In this section we collect the results for, in turn, the 
fluctuation corrections to the sound attenuation and then 
to the sound velocity. 



A. Sound attenuation 

The odd powers of the w-cxpansion of the analytically 
continued results for the diagrams of Fig.Elgfve the fluc- 
tuation corrections to the imaginary part Im (ll R (uj)) of 
the retarded phonon self-energy. One can then use this 
result to obtain the sound attenuation from the imagi- 
nary part of the pole of the propagator, as in Eqs. © 
and 10. 

All diagrams contribute to the polarization propagator 
to leading order in u. The different contributions to the 
sound attenuation coefficient arc then extracted from the 
imaginary part of the retarded fluctuation polarization 
operator. In the most general form, they can be written 



a^(T,u) = ^^(TrJ/Wfcr.T,), (28) 
e F v s 

where (3 denotes the particular channel (DOS, rMT, aMT 
or AL). In Eq. /^(e, r,7^) means the function of 

temperature which comes from the integrals in Eqs. (|20|l . 
(12211 . 1231) . and (|27|l . and which also depends on the mate- 
rial properties. The coefficient k^] (Tt) comes from the 
integration of Green functions and impurity blocks in the 
expression of the type l|16[) . It has weaker temperature 
dependence at e — ► than f^'(e,r, 70), and it basically 
shows how the impurity concentration affects the result. 

In the following, it will be more convenient to use the 
dimensionless parameter t)(Tt) =r]/ (tuf) 2 . We list here 
the functions which also appear in the conductivity 



calculation, 122] and the coefficients k 



(JJ) 



For the reader's 



convenience, we also give the limiting expressions of 
which depend on the relation between e, r, and 7^, where 
applicable. 
DOS: 



/( DOS )(e,r) = ln 



In r , e <C r, 
lne , r <C e. 



ados),, _ £ (|) - W (| + M 

K att ( X ) — 



(29) 



8ir 2 xrj(x) 
28C(3) 



3^2 



, X < 1, 



(30) 



7C(3) ' 



x > 1. 



Regular part of MT: 



(31) 




r < e. 



^f T) (x) 



n 3 x - 2nxip' (\ + j^) + xjj" (\\ 

ia>-K 2 xf){x) 

14C(3) 

^3 , sc<l, 



(32) 



2x, 

Anomalous part of MT: 



x > 1. 



/(° MT )(e,r,7 ) 



1 



1 



70 



2e ' 

1 

270 



+ O + V /e ( e + r ), 

e < 70 >C r, 

70 < e < r, 

70 < r < e, 
e < r < 70, 



(33) 



, 87^' ^«min{70,e}. 



(aMT) _ _7T 
^att — ^ • 



(34) 
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AL: 



VI. DISCUSSION 



(AL,a) 



1 

r ' 



2(e + r) 



3r_ 

8e 2 ' 



e < r, 
r < e. 



(35) 



(36) 



The fluctuation terms scale like w 2 just as in the corre- 
sponding normal metal case without fluctuation correc- 
tion. 



B. Sound velocity 

The renormalization of the phonon frequency w(k) ob- 
tained from Eq. IjlUfl provides the fluctuation corrections 
to the sound velocity when we use oj(k) = v s k. It should 
be emphasized that there is no anomalous MT term in 
leading (zcroth) order of lj 1 while to higher order of to all 
types of temperature dependences appear. To leading or- 
der uj° the non-zero terms give, for Auj/uj = (lu — loq) /ujq, 
the expressions in a form: 



, lo) Avs (T, u) 
uj v s 

£F 



(37) 



For DOS, rMT and aMT terms, the temperature func- 
tions / (/3) ( e ; r ;7</>) are S iven b y Eq s - CHJl and (J33J), 
respectively. The temperature function f^ AL,v '(e,r) for 
the AL diagram is given by the Eq. (|41|) . which is differ- 
ent from that for attenuation (|35() . We then obtain for 
the coefficients k^: 



1 



DOS: 

(DOS) _ 
K v — 1, 

Regular part of MT: 

(rMT) 

«„ - 2 
Anomalous part of MT: 

4° MT) = 0, 

AL: 

K (AL) f (AL, v){ejr) 



(38) 
(39) 

(40) 
(41) 



1 , e < r, 
3r 



Comparison of the normal state results without fluctu- 
ation corrections found in Sec. [H] with Eqs. I|28[) and l|37|) 
above show that the fluctuation corrections arc smaller 
than the normal state contribution by a factor of T/ep. 
For layered organics, the Fermi surface parameters of Ref. 
[25j lead to T/sf ~ 1CP 2 . The temperature functions 
f(P> increase this ratio at T — > T c , thus making the fluc- 
tuation corrections experimentally measurable. 

It can be seen from the results of previous section that 
it is the k^'(Tt) coefficients that determine the sign 
of the fluctuation sound attenuation and sound veloc- 
ity. The coefficients for sound attenuation fi4f](TV) in 
Eq. I|28|) have weaker temperature dependence than that 
of the functions f^. One can see that these coefficients 
in all fluctuation diagrams are finite as (Tt) — > 0, and 
that they are of the same order of magnitude in this 
limit. In the opposite (clean) limit (Tt) 3> 1, these coef- 
ficients increase as a power law of (Tt). We recall that 
there is an analogous formal divergence of the DOS and 
MT coefficients [22| in the fluctuation conductivity but in 
that case they cancel each other in the limit (Tt) ^> 1. 
In our case, the formal extension of our results to the 
limit (Tt) ^> would be incorrect because the local 

approximation for the fluctuation propagator (fill) and 
the impurity vertex (|13f) have a natural limit of validity 
(Tt) <C (the so-called local clean case). |2|| 



I 4e 



r < e. 



In contrast to the coefficients for sound attenuation 
k~m(Tt), the coefficients in the case of sound velocity, 
Kv in Eq. (|37fl . do not depend on Tt, thus simplify- 
ing the estimation of the sound velocity corrections in 
particular materials. 

The overall magnitude of the fluctuation correction is 
given mainly by the temperature functions These 
functions depend on the anisotropy parameter r and the 
phase-breaking time r^, which are defined by the intrin- 
sic material properties. The temperature function corre- 
sponding to a given type of diagram is identical for sound 
attenuation and velocity, except for the AL contribution 
which has temperature dependences that differ. In gen- 
eral, all the fluctuation terms remain finite as T — > T c . 
The asymptotics of the corresponding temperature func- 
tions have been given in the preceding section and they 
were also discussed elsewhere. 

The functions f^' also describe the crossover between 
two and three dimensions, where by three-dimensions 
we mean that the superconducting correlation length 
is much larger than the interlayer spacing. The two- 
dimensional behavior corresponds toe > r and the three- 
dimensional behavior to the opposite limit r 3> e. The 
crossover occurs at e ~ r, namely when the supcrcon- 
ducting correlation length £j_(0) \e\ in the perpen- 
dicular direction becomes of the order of the interplanc 
spacing. 

In order to compare the contributions from different 
diagrams, we consider a range of parameters which is 
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0.000 0.002 0.004 0.006 0.008 0.010 alistic 7^ = 0.1 and at various values of anisotropy parameter 

r = 1CP 1 . . . 1CP 4 , b) at most realistic for the k-(ET)2 family 
c c compounds r — 10~ 4 and at various 7^ = 10 _1 . . . 10 -6 . 



FIG. 4: The temperature dependence of f^ r ' at various 
values of anisotropy parameter r = 1CP 1 . . . 10 -4 . 



realistic for layered organic superconductors. The esti- 
mation for the anisotropy t\\/t± « 4000, where ti is the 
intralaycr transfer integral, is in agreement with the mag- 
netoresistance data for k-(ET) 2 Cu(NCS) 2 . H3 This gives 
an anisotropy parameter of the order of r > 10 -4 . 

Concerning the scattering time, its estimate |28| from 
Shubnikov-de Haas experiments in k-(ET)2 materials 
gives t w 3 ps, which leads to Tt ~ 1 at T > T c . 

The phase-breaking time can be estimated from the 
experimental data on fluctuation contribution to trans- 
port coefficients using the appropriate expression for the 
MT diagram. The analysis of the magnetoresistance data 
on YBCO in Ref. H3, which included AL and MT con- 
tributions, resulted in values ~ 1/T. Later on, the 
authors of Ref. (2^ claimed that neglecting the orbital 
pair breaking effects while leaving the Zceman contribu- 



tion to the MT diagram in Ref. can be incorrect. 
Hence, they extended the range of the realistic values of 
to T0 < 10/T with ~ 100/T being still acceptable 
for the analysis. l22l A similar estimate was chosen by the 
authors of Ref. [T(| to predict the magnitude of the fluc- 
tuation corrections to the nuclear spin-lattice relaxation 
rate and the NMR Knight shift in high-T c cuprates. To 
the best of our knowledge, no relevant analysis in layered 
organic materials have been reported to date. However, 
it is known that the phase-breaking time should satisfy 
the relation > r. Without loss of generality then, we 
assume that pair breaking by inelastic scattering is weak 
enough and we consider 7^ ~ 10 -3 ...10 -2 . 

The aMT contribution vanishes to leading order, 
uj°, for sound velocity. With the above estimation 
for the anisotropy parameter, one can see that the 
AL part as well as the regular MT part give the 
least important non-vanishing contributions to sound 
velocity (order Indeed, lim 6 _^o f^ rMT H £ i r ) ~ 

\im £ ^of {AL ' v) (e,r) = 1 at any r, while for the DOS con- 



10 



10000 



1000 




0.004 0.006 



0.010 



b) 1.0 

0.8 
0.6 
0.4 
0.2 
0.0 



-r=10 
- r=10" 3 
r=10" 2 
r=10" 1 



0.000 0.002 0.004 0.006 0.008 0.010 
(T-T c )/T c 



FIG. 6: The temperature dependence of AL contributions 
a) to the sound attenuation f' L ' a ' , b) to the sound veloc- 
ity f( AL ' v ^ a t various values of anisotropy parameter r = 

io _1 . . . i(n 4 . 



tribution, lim e _ f {DOS) (e, r) = (1/2) ln(4/r) w 5 at 
r ~ 10" 4 . 

For sound attenuation the contribution of the AL term 
can be significant near T c , by contrast to its contribution 
to the sound velocity, since hm s — > of^ AL ' a ^(e,r) = l/ r - 
However, it decays much faster than all other fluctuation 
terms above T c : f AL > a )( e = 0.002)//( Ai > Q )(e = 0.01) ~ 
23 with lim^/e^o f AL ^\e) ~ 1/e 2 . At the same time, at 
70 ~ 10~ 3 ...10~ 2 and r ~ 10 -4 , we see that the anoma- 
lous MT temperature function is large, f( aMT >(e = 0) w 
500 at 7 - 10~ 3 and f aMT ){ e = 0) w 50 at 1<j} ~ 10~ 2 . 
Thus, at some 7^, the aMT can be much larger than all 
other contributions. However, while it is quite large at 
e = 0, the aMT part decreases with increasing T faster 
than the DOS contribution. Indeed, at r = 10~ 4 , the ra- 
tio f( DOS \e = 0.002)//( £los )(e = 0.01) ~ 1.3 shows 
that the DOS contribution has not changed much in 
that temperature range while the corresponding ratio for 
aMT is f( aMT )( e = 0.002)//O MT )(e = 0.01) ~ 5 for the 
same r and in the relevant range of 7^. In other words, 



the logarithmic asymptotics l|29|) of the DOS diagram at 
r/e — > wins over the power-law dependencies of the 
aMT term 

One can see that next order terms in the frequency 
expansion for the sound velocity and sound attenuation 
would contain the factor (ui/T) 2n , where n = 1,2,.... 
The typical ultrasound frequency is of the order of u> ~ 
100 MHz. At T > T c ~ 10 K, the relevant frequency to 
temperature ratio is thus of the order (u/T) 2 ~ 10~ 7 . 
One can check, for example, that in order to compensate 
for such a small factor, the parameter 7^ in the anoma- 
lous MT term at r ~ 10 -4 should be as small as 10 -10 , 
corresponding to unreasonably large ~ 10 8 r. We con- 
clude that, in our model, the sound velocity should be in- 
creased by the fluctuation corrections given by the DOS 
term. 

Finally we would like to comment on the applicability 
of our model with regard to the character of interlayer 
electron transport in organic conductors, a subject which 
is still hotly debated. The interlayer tunneling of elec- 
trons can be coherent or incoherent (see Refs. [3CX l3ll | and 
references therein). Following the terminology of Ref. 
|3lj . we should distinguish between two limiting cases 
of incoherent interlayer transport, strongly and weakly 
incoherent. First, in the strongly incoherent case the in- 
tralayer electron momentum is not conserved in the tun- 
neling processes between adjacent layers because the tun- 
neling can be accompanied by strong elastic or inelastic 
processes. In this case the motion of electrons between 
the layers is diffusive and the electron band states in the 
z direction and the corresponding Fermi surface cannot 
be defined. Our model is invalid in that limit. It is valid, 
however, in the second case, namely for weakly incoher- 
ent tunneling, that occurs when rt± <C 1. In that case the 
intralayer electron momentum is conserved in the tunnel- 
ing process and the electron wave function in adjacent 
layers has some overlap, but there are many in-plane col- 
lisions between tunneling events. Our approach is also 
valid obviously in the coherent case where rt^ 1. 

It is clear from our results that we can also consider 
two other limiting cases, namely the clean limit (TV 3> 1) 
and the dirty limit (TV 1). The latter always corre- 
sponds to weakly incoherent tunneling while both coher- 
ent and weakly incoherent tunneling can occur in the 
former. This can be explained as follows. We have as- 
sumed throughout that r <S 1. That parameter takes the 
following limiting forms: 



4£(0) 



T 



T 2 



, Tt < 1 



Tr> 1 



(42) 



where at T = T c , £j_(0) is the size of the Cooper pair 
perpendicular to the plane. The conditions r«l , and 
Tt <C 1, corresponding to dirty case, are consistent only 
with the weakly incoherent tunneling since they lead to 
(Tt±_) 2 < Tt < 1. On the other hand, for Tt > 1 the 
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condition r <C 1 is consistent with both coherent and 
weakly incoherent tunneling. This explains why in the 
preceeding discussion we distinguished only between the 
clean and dirty limits. Note that the crossover from clean 
to dirty limit occurs at vf/T ~ VpT, namely when the 
elastic mean free path becomes smaller than the thermal 
dc Broglie wavelength. 

In summary, the superconducting fluctuation correc- 
tions to the sound attenuation, in a realistic range of pa- 
rameters, are given by the sum of the DOS, anomalous 
MT and AL terms, which decrease the normal state at- 
tenuation. In contrast, to leading order in to, the correc- 
tions to the sound velocity are given by DOS diagrams, 
because the expansion for the anomalous MT diagram 
begins at order to while the AL diagram is small to this 
order. 



VII. CONCLUSION 



the sound velocity as we approach T c . The anomalous 
contribution of the Maki-Thompson diagram is propor- 
tional to (oj/T) 2 and generally can be neglected unless 
phase breaking is extremely low. Further details of our 
calculations may be found in Ref . ■ 

It should be noted that our results on sound attenu- 
ation are in qualitative agreement with the preliminary 
experimental data. Q These experiments are currently in 
progress. The detailed analysis including the comparison 
with the present theory will be published elsewhere. 

To make actual comparisons with experiment, one may 
need to consider several other physical effects. For exam- 
ple, strongly incoherent interplane tunneling |3lj | and d- 
wave superconductivity [3j| may have to be considered. 
In addition, one may need to include in the fluctuation 
propagator Eq. (|11J) mode-coupling terms that will allow 
to go beyond mean-field theory in the description of the 
crossover from two to three dimensions. 



We calculated, in the hydrodynamic limit lot -C 1 and 
H<1, the effect of superconducting fluctuations on lon- 
gitudinal sound propagating perpendicular to the layers 
of quasi two-dimensional systems above T c . This corre- 
sponds to the experimentally realizable situation. The 
detailed results for the sound attenuation coefficient a 
and the sound velocity renormalization Av s /v s are sum- 
marized in section ^ In short, there are four types 
of temperature dependences coming respectively from a) 
the density of states diagrams, b) the regular contribu- 
tion of the Maki-Thompson diagram, c) the anomalous 
contribution of the Maki-Thompson diagram and d) from 
AL diagram. These contributions have different signs 
and their magnitudes strongly depend on material prop- 
erties such as anisotropy and phase-breaking processes. 
For longitudinal sound propagating perpendicular to the 
layers, the open Fermi surface condition leads to a sup- 
pression of the divergent contributions to leading order, 
so that fluctuation corrections to sound attenuation are 
expected to generally be much smaller than the fluctua- 
tion corrections to conductivity. 

The leading temperature dependences with respect to 
T c of the fluctuation corrections to the sound attenuation 
have the same to 2 scaling with frequency as in the impure 
(dirty) metal without fluctuation correction. To this or- 
der in to, all four types of temperature contributions are 
present. As we move close to T c , the fluctuation contri- 
bution decreases the sound attenuation below its normal 
state level, at least for realistic values of the material 
parameters. 

The leading fluctuation correction to the sound veloc- 
ity does not depend on to, just like the normal state term. 
To this order, there are only three types of temperature 
dependences coming from the density of states diagram, 
AL diagram and from the regular contribution of the 
Maki-Thompson diagram. The overall fluctuation correc- 
tion in realistic material situation should be dominated 
by the density of states contribution and should increase 
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APPENDIX A: ON THE ANALYTICAL 
CONTINUATION OF THE DIGAMMA 
FUNCTION OF COMPLEX ARGUMENT 

We shall discuss an important mathematical issue oc- 
curring in the evaluation of the anomalous MT contribu- 
tion i|18|) and its further analytic continuation. Namely, 
let us look at the typical sum over the fermionic Matsub- 
ara frequencies: 

Z=l 1 

r£-. (Ai) 

71=0 £ ™ 

The expression (|A1|) appears if one decomposes into par- 
tial fractions the initial sum over e n in the anomalous 
part of the Maki-Thompson diagram and neglects the 
small terms ~ Dq 2 in them. One can easily evaluate 
the sum (lAlf) using the well-known representation of the 
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digamma function: 

^ = ^K^)+21n2 + 7 }, (A2) 

n— v \ / j 

where 7 = 0.577216 is the Euler constant. Here v 
comes from the external bosonic Matsubara frequency 
u> v = 2itTi>, with v an integer. At the same time, it is 
easy to see that the straightforward change of variable 
n = —n' — 1 transforms the initial expression IjAljl in a 
way that one obtains the alternate result: 

v-l -1 

r Ef = - T E^- ( A3 ) 

7i—0 n' — — v " 

= i{H^ iy ) +21n2 + 7 }- 

As long as we deal with the Matsubara frequencies uj u = 
2irTv these two results Eq. (|A2|) and Eq. <| A3|> are equiv- 
alent because of the relation 

^(l~ z J =^(l +z ) +7rtan(7rz), (A4) 

and the fact that tan(7rz) = when z = v is an integer. 
Therefore, at this stage we are free to choose either Eq. 
(|A2() or Eq. I|A3[) as the result of the summation. The 
distinction appears when one tries to do the analytic con- 
tinuation to real frequencies iLo u — > Lu + irj. After analytic 



continuation, Eqs. (|A2|) and (|A3|) . which we obtained 
from the same diagram, are not equal to each other any- 
more. Indeed, for non-integer v, ip fg — 2) =Jp(h z ) 
is not true anymore, as one can see from Eq. (|A4|) . 

It is crucial to choose either Eq. (|A2|) or Eq. (|A3|) 

as the result of the summation before we do the ana- 
lytic continuation. The choice is determined by whether 
we want the advanced or the retarded response func- 
tion. Indeed, consider Eqs. 1A2|I and l|A3() as general 
response functions \ M i^v) = x{ z )\z=iuj v ^ ip (§ ± v) on 
the Matsubara frequency domain uj u . If we then ana- 
lytically continue the argument z = iuj v — ► cu + iij of 
ip (1/2 — iz/(2ttT)), the poles of resulting function arc 
all located below the real axis of the complex plane thus 
giving us the retarded function On the other 

hand, the analytic continuation z = iui v — ► uj — ir] of the 
argument of ip (1/2 + iz/(2nT)) gives us the advanced 
response because all its poles are located in the upper 
half-plane. 

We thus conclude that if the digamma function appears 
in our calculations, we are free to use either Eqs. (|A2|) 
or i A3|) because they are equivalent, as long as we deal 
with Matsubara frequencies. However, before doing the 
analytic continuation, we have to decide whether we need 
the advanced or the retarded function. If we want the 
retarded function, we use the rule given by Eq. (jJ2j and 
deal only with digamma functions of the type ip(l/2 + 
v). Otherwise, use ip(l/2 — v) to obtain the advanced 
response. 
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